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ABSTRACT 


The experimental study of the interaction of atomic oxygen with organic polymer 
films under LEO conditions has been hampered by the inability to conduct detailed 
experiments in situ. As a result, studies of the mechanism of oxygen atom reactions have 
relied on laboratory O-atom sources that do not fully reproduce the orbital environment. 
For example, it is well established that only ground electronic state oxygen atoms are 
present at LEO, yet most ground-based oxygen sources are known to produce singlet 
oxygen atoms and molecules and ions in addition to 0( 3 P). Engineers should not rely on 
such facilities unless it can be demonstrated either that these different oxygen species are 
inert or that they react in the same fashion as ground state atoms. 

Ab initio quantum chemical calculations carried out at NASA Ames and elsewhere 
have been aimed at elucidating the biradical intermediates formed during the electrophilic 
addition of ground and excited-state oxygen atoms to carbon-carbon double bonds in 
small olefins and aromatic molecules. These biradicals are critical intermediates in any 
possible insertion, addition and elimination reaction mechanisms. Through these calcula- 
tions, we will be able to comment on the relative importance of these pathways for 0( 3 P) 
and 0( 1 D) reactions. The reactions of oxygen atoms with ethylene and benzene will be 
used to illustrate the important features of the mechanisms of atomic oxygen reaction 
with unsaturated organic compounds and polymeric materials. 
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I. INTRODUCTION 


The kinetics and mechanisms of the chemical reactions of oxygen atoms with organic 
molecules have been the subject of intensive study for more than thirty years. Most 
of these investigations have been carried out in the gas phase as the major impetus 
for this work has come from the combustion and atmospheric photochemistry research 
communities. Cvetanovic and Singleton(l) have recently reviewed this work for the case 
of O-atom olefin reactions. Companion studies have been carried out by Klein and 
Scherr(2) in low temperature glasses. These studies have been fraught with difficulties, 
as it is not an easy task to generated a clean, controlled and calibrated source of oxygen 
atoms. More recently, molecular beam studies of 0( 3 P) reactions with organic species 
under single-collision conditions have provided a detailed look at the primary processes in 
these complicated reactions (3-6). A comprehensive reaction mechanism must be capable 
of describing these single-collision, gas phase and condensed phase experiments. 

The chemistry of atomic oxygen is made especially interesting by that species’ triplet 
ground electronic state and by the existence of a low-lying excited singlet state(7). The 
following general conclusions can be drawn from this large body of work (1,8, 9): 0( 3 P) 
abstracts hydrogen from saturated organic molecules and adds to carbon- 
carbon double bonds resulting in the formation of carbonyl compounds and 
epoxides, while 0( 1 D) tends to form alcohols by inserting into C-H bonds. 
These initial products undergo secondary reactions which lead to fragmentation of the 
reactant molecules. In most cases, the branching ratios for secondary fragmentation prod- 
ucts are statistical, meaning they are governed by the relative thermodynamic stabilities 
of the energetically permissible products (or by the relative ease of bond rupture). 

Renewed interest in the kinetics and mechanisms of atomic oxygen reactions has been 
generated during the early flights of the space shuttle by observations that most organic 
polymeric materials undergo noticeable degradation during exposure to the low Earth 
orbit (LEO) environment(10,ll). The most severe damage was observed on surfaces 
facing the direction of motion in orbit. In this orientation, both the flux and the kinetic 
energy of the atmospheric species striking the surface are high (« 10 15 cm -2 s -1 and 5-10 
eV). Ground state atomic oxygen is the principal constituent of the residual atmosphere 
in the 150-250 km altitude range and it seems likely that this species would react readily 
with most organic materials. Dedicated flight experiments on shuttle flights STS-5 and 
STS-8 have confirmed these early observations (10), indicating that the degradation is 
widespread and independent of surface temperature and solar UV/VUV irradiation. 

Several laboratory facilities (12) have been constructed in an attempt to simulate 
the LEO environment for kinetic and mechanistic studies and for screening materials 
proposed for use on the space telescope and space station projects. Unfortunately none 
of these facilities is capable of completely achieving the desired conditions: a beam of 5 
eV 0( 3 P) atoms of the appropriate flux and without contamination from O 2 , 0" r , 0( J D) 
or other reactive species. 
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As one purpose of this workshop is to assess the current status of laboratory simula- 
tion of the LEO environment, it is appropriate to examine the differences in the reaction 
mechanisms of 0( 3 P) and O^D) with organic compounds. It should be pointed out 
that very little is known about O atom reactions with polymers (13) and there is no cer- 
tainty that the same mechanisms apply to the reactions of small molecules and polymers. 
However, it is likely that the initial steps in the reaction mechanisms are the same in 
both cases. In the present study we consider the reactions of atomic oxygen atoms with 
unsaturated organic molecules from a theoretical point of view. We describe a previous 
computational study of the addition of oxygen to ethylene (14) and present the results of 
new calculations for the reaction of 0( 3 P) + benzene. These systems are chosen because 
of the availability of experimental data and because they have features in common with 
organic polymers being considered for LEO applications such as KAPTON (a polyimide 
film) and structural materials such as PEEK and polysulfone. Some of the questions we 
will attempt to answer are: “under LEO conditions will 0( 3 P) and 0( J D) react 
differently?”, “are the kinetics affected by the 5 eV O atom energy?”, and 
“how do the degradation rates scale with Q-atom flux?”. 


II. COMPUTATIONAL DETAILS 


In recent years, computational chemistry groups at NASA Ames Research Center 
and elsewhere have been able to capitalize on the availability of supercomputers such as 
the Cray XMP and have been performing accurate “first principles” calculations of the 
spectroscopic, kinetic and thermodynamic properties of small molecules. These data arise 
from solutions of the Schrodinger equation which describes the energy and distribution 
{wave function ) of a collection of electrons in the field created by the motionless nuclei 
of the atoms in a molecule(15). By relaxing the locations of the atoms in a series of 
calculations, molecular equilibrium geometries and saddle points for chemical reactions 
can be determined. Collisional and spectroscopic properties are obtained by subsequent 
calculations in which the nuclei are allowed to move, subject to the forces generated by 
the gradient of the electronic energies. 

The calculations reported herein utilize several “standard” aspects of ab initio quan- 
tum chemistry codes(15). The molecular wavefunction is expanded in a basis of Gaussian 
functions which represent atomic orbitals and, at the Hartree-Fock or Self-Consistent 
Field (SCF) level of calculation, each electron experiences only the average field of the 
remaining electrons in the molecule. SCF calculations give only an approximate picture 
of the electron density distribution in a molecule. A more accurate description is obtained 
by allowing each electron to feel the instantaneous influence of all the other electrons. 
Various computational schemes have been developed to include these effects, two exam- 
ples of which are the multiconfiguration self-consistent field (MCSCF) and configuration 
interaction (Cl) methods. These effects are generally of critical importance in the study 
of chemical bond breaking and formation. However, experience has shown that SCF 
calculations are quite adequate for the determination of geometries for stable molecules 
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and reaction intermediates. 


In the present study we have used all of the above mentioned methods. The calcu- 
lations were carried out using Dunning’s double-^- atomic orbital basis sets in which 9s 
and 5p Gaussian functions are contracted to make 4s and 2p orbitals for each carbon and 
oxygen atom and 5s Gaussians are contracted to 3s orbitals for each hydrogen(16). The 
molecular geometries were determined by minimizing the norm of the analytic gradient 
of the SCF energy using a computer code GRADSCF(17). At these stationary points, 
the analytic second derivatives of the energy were evaluated to determine whether the 
computed geometries correspond to local minima or to saddle points(18). At the SCF 
geometries, a variation of MCSCF calculation known as the Complete Active Space SCF 
method (CASSCF) was employed(l9). In these calculations, the oxygen 2p electrons and 
orbitals and the benzene 7r electrons and orbitals were included in the “active space” 
and the interactions between all possible arrangements of those 10 electrons in the 9 or- 
bitals were considered. Finally, contracted Cl calculations(20) were performed using the 
CASSCF orbitals as reference vectors and including single and double excitations from 
a set of reference configurations comprised of all configurations making important con- 
tributions to the CASSCF wavefunction at any of the molecular geometries considered. 
For O -f benzene in Cs symmetry 3 A\ 3 A", *A' and 1 A" states are possible. The sizes 
of the CASSCF and CCI wavefunctions generated for each of these states are given in 
Table I. 


III. Q( 3 P) + ETHYLENE 


The reaction of ground state oxygen atoms with ethylene has been the subject of 
extensive experimental(l,3-6,8,21-23) and theoretical(l4) study. The gas phase kinetics 
studies demonstrated that 0( 3 P) readily adds to the carbon-carbon double bond, re- 
sulting in the formation of a transient biradical species which undergoes triplet-singlet 
interconversion and rearrangement leading to epoxide and aldehyde products. However, 
molecular beam experiments under single-collision conditions found rapid hydrogen atom 
elimination follows the biradical formation almost exclusively (4-6,23). The elimination 
pathway is favored under these conditions because it does not require a change in electron 
.spin. Carbon-carbon bond rupture to form fragmentation products (either formaldehyde 
and methylene or, after a 1,2-hydrogen shift, formyl and methyl radicals) is possible, 
but does not seem to occur. Similarly, the rearrangement to form triplet acetaldehyde is 
not observed, probably because the hydrogen shift involves surmounting a higher energy 
barrier than does hydrogen elimination. Figure 1 is a schematic energy diagram (taken 
from ref. 8) which describes these processes. 

One of the important results of the theoretical study is that the singlet and triplet 
biradicals are nearly degenerate and intersystem crossing should occur readily under mul- 
tiple collision conditions (e. g., bulk gas or condensed phase). Once the system has been 
converted to the singlet manifold, rearrangement to form ethylene oxide or acetaldehyde 
can compete with the still-present elimination channel. The addition products are ex- 
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tremely stable compared to the biradical (« 400 kJ/mol for acetaldehyde), but collisional 
stabilization is required to prevent their unimolecular fragmentation. 

The reactions of 0( X D) can proceed by the same biradical-addition pathway as in 
the triplet case, or by direct insertion into the C-H bond to form (C 2 H 3 )OH. The latter 
channel generally dominates the kinetics, however, because it involves much stronger 
attractive forces between the reactants. The resulting alcohol will undergo fragmentation 
unless it is collisionally stablized. 

Hydrogen abstraction is also possible in the reaction of atomic oxygen with unsat- 
urated compounds, but its occurrence is negligible at room temperature(24). However, 
at elevated temperatures, it is expected to become increasingly more important. It is 
possible that it can even compete with the insertion process in the O^D) chemistry (9). 

The abstraction and elimination pathways lead to the formation of free radical prod- 
ucts that will certainly undergo further bimolecular reaction. The unimolecular fragmen- 
tation products will also be highly reactive. All of these fragmentation processes will 
tend to produce statistical product distributions, with the weakest bonds being the ones 
most frequently broken. The major difference between ultimate abstraction, elimination, 
addition and insertion products will most likely be in the amount and nature of the 
oxygen content. Hydrocarbons and water would be formed from abstraction, carbonyl 
compounds would be formed from elimination and addition, and hydroxyl compounds 
from insertion. 


IV. Q( 3 P) + BENZENE 


The reactions of 0( 3 P) with benzene and toluene have also been the subject of con- 
siderable study(8, 25-29). The measured activation energies(8,25,26), 15-20 kJ/mol, are 
intermediate between the values found for addition to olefins and hydrogen abstraction. 
For temperatures as high as 950 K(26), no curvature was observed in the Arrhenius plots 
of the rate constant data and no hydrogen/deuterium isotope effect was seen. Molec- 
ular beam experiments (28,29) have clearly demonstrated that elimination products are 
dominant under single-collision conditions at collision energies between 10.5 and 41.0 
kJ/mol. Curiously, however, for the case of 0( 3 P) + toluene, Lee and coworkers(29) 
mainly find methyl elimination products (i.e., phenoxy and methyl radicals). It seems 
unlikely that a localized biradical adduct should form at that site, exclusively. However, 
the methyl-phenvl bond energy is about 38 kJ/mol less than the H-phenyl bond energy. 
Perhaps facile 1,2-oxygen shifts take place until the O-atom migrates around the ring 
to the methyl carbon, permitting cleavage of the weaker methyl-phenyl bond to occur. 
This 1,2-oxygen shift would involve epoxide-like intermediate structures. As the ab ini- 
tio quantum chemical study of the ethylene-oxygen adduct(14) provided insight into the 
mechanism of the 0( 3 P) + ethylene reaction, we have undertaken a similar study for O 
+ benzene(30). 
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The initial calculations for the 0( 3 P, 1 D) + benzene system involved determination 
of the equilibrium geometries for the reactant (benzene), biradical adducts, the addition 
product (epoxide) and insertion product (phenol). SCF gradient calculations were per- 
formed for the closed shell benzene, epoxide and phenol molecules, for the triplet and 
singlet open-shell biradical adducts and for the triplet epoxide intermediate. The geom- 
etry of the lowest energy biradical ( 3 A') is shown in figure 2. It has one radical electron 
in what is mainly a benzene 7r orbital (16a') which is localized on the ortho and para 
carbons. The other unpaired electron (17a') is localized on the O-atom and lies in the 
O-C-H plane perpendicular to the C-G bond. The Q-C-H bond angle is 103° with the 
O and H atoms evenly distributed above and below the benzene ring plane. The carbon 
atom involved in the adduct bond is nearly tetrahedral and the C-0 and C-C bonds 
are rather elongated. The individual CASSCF orbital densities are plotted in figure 3 
for some of the triplet biradical ( 3 A') valence orbitals. It can be readily seen that the 
aromaticity is lost and that one of the unpaired biradical electrons occupies a benzene ■n 
orbital while the other is localized on the oxygen atom. 

The resulting energies from the SCF, CASSCF and CCI calculations are summarized 
in Table II. It can be seen that significant differences exist between the experimental en- 
ergy estimates and the results of the large (MR)-CCI calculations. The approximation 
that only the benzene n electrons and the O-atom 2p electrons have been correlated is 
more severe for phenol with its O-H bond than for the biradical. Also a larger basis 
set would have to be used for more accurate calculations. As was observed for O + 
ethylene(14), the singlet-triplet energy splittings in the biradicals are quite small. This 
suggests that interconversion would occur fairly rapidly. The SCF calculations for the 3 A' 
and 3 A" biradicals and the 3 A" epoxide included determination of the harmonic vibra- 
tional force fields and normal mode frequencies. Examination of those results indicated 
that the structures were all true local minima(18), implying the existence of barriers 
to the interconversion of biradicals to epoxides or phenols. The locations and heights 
of these barriers have not been determined in the present study. However, the triplet 
epoxide is considerably higher in energy than the biradical which makes the possibility 
of 1,2-oxygen shifts seem unlikely. 

The similarity between the energetics for O-atom reactions with ethylene and ben- 
zene (or toluene) indicate that the kinetics should also be similar. One expects elimina- 
tion to dominate under single-collision conditions with addition being most important for 
bulk gas or condensed phase experiments. At higher energies, hydrogen abstraction would 
also be expected to play an increasingly more important role. The difference in products 
between 0( 3 P) and 0( J D) chemistry under multiple collision conditions is not readily 
apparent, as the primary addition and insertion products would all undergo secondary 
fragmentation. 


V. IMPLICATIONS FOR POLYMER REACTIONS 


It is tempting to extrapolate the reaction mechanisms described above to the case of 
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O + polymer reactions under thermal and high energy conditions. One expects the initial 
reaction steps to involve similar abstraction, addition, elimination or insertion processes. 
However, a polymer can be thought of as the ultimate condensed phase molecule and 
rapid stabilization and internal conversion should be possible. In general, fragmentation 
products should be formed based on their thermodynamic stability and the relative weak- 
ness of the bonds being broken. For the case of polystyrene, one might expect to observe 
a loss of hydrogen atoms (due to abstraction) and phenyl groups (due to elimination and 
fragmentation) from the chain backbone. The resulting alkyl radicals would probably 
behave similarly to the intermediate formed in the reaction of 0( 3 P) with polyethylene. 

If translationally hot oxygen atoms were present, small energy barriers would not 
impede reaction and the abstraction process would probably be enhanced compared to 
addition or elimination. This change in reaction mechanism means that rate constants 
cannot readily be extrapolated from the thermal to 5 eV O-atom case. Furthermore, 
the products for 0( X D) reactions would be similar to those for ground state oxygen at 
higher energies (the reactions are less discriminating) even though they are quite different 
at lower energies. Also, it is possible that the magnitudes of the rate constants for the 
0( 3 P) and 0( ! D) reactions would be similar at high energies if both processes are limited 
mainly by the collision rate. 

Two other critical factors depend on the condensed-phase nature of the polymer 
films. If the oxygen atoms are accomodated (t. t., thermalized) on the solid surface 
before reacting, the 5 eV translational energy will not be an important factor and the 
nuances of the thermal reaction mechanisms will be important. Secondly, under the 
high vacuum conditions encountered in both LEO and molecular beam experiments, 
volatilization of moderate molecular weight oligomers will occur. If the O + polymer 
reaction results in a downward shift in the molecular weight distribution of the polymer, 
enhanced vaporization will take place which will result in a significant contribution to the 
measured weight loss in vacuum, but not at higher pressures. These concerns, which are 
beyond the scope of the present theoretical study, need to be addressed if a comprehensive 
description of atomic oxygen-polymer reactions is to be developed. 
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Table I. Size of the Calculations 


Symmetry 0, 

CASSCF 6 

•^config. 

(MR)-CCI 

N ?ef. 

•^config. 

*A' 

2688 

6 

368,771 

1 A" 

2604 

12 

792,295 

3 A' 

3756 

9 

1,143,225 

3 A" 

3804 

13 

2,225,971 

3 pjtd 

3800 

8 

573,296 


° Geometries restricted to C s symmetry (O and two para C atoms lie in the symmetry 
plane) . 

6 9 orbitals (6 a' and 3 a”) and 10 electrons included in the active space (includes O 
2p shell and benzene n orbitals) . 

c Number of reference configurations used in the CCI. For each symmetry and multi- 
plicity this includes all CASSCF configurations with coefficients greater than 0.075 
for any geometry. 

d Epoxide calculations (symmetry plane bisects 2 C-C bonds) based on 5 a' and 4 a" 
■orbitals in the CASSCF active space. 


83 



Table II. Computed Energies 0 for CeHe + O ( 3 P and *D) 


Symmetry 

SCF 6 

CASSCF C 

(MR)-CCI 

EXPT. d 

Asymptote 





3 A" 

0 

0 

0 

0 

r A " 

215.9 



188.7 

Biradical 





3 A' 

151.0 

41.8 

9.6 

(-63.) 

3 A" 


54.0 

53.6 

J A' 

156.9 

51.5 

22.2 


1 A" 


59.4 

57.7 


Phenol 





J A' (So) 

-106.3 

-228.4 

-280.7 

-428.4 

3 A' (TO 

303.3 

156.1 

100.8 

-87.9 

! A' (Sj) 




14.2 

3 a" 


274.1 

196.2 


l A" 


266.1 

202.5 


Epoxide 





*A' (So) 

138.5 




T— 1 

CO 

306.7 

166.1 

144.3 



° Energies in kJ/mol. 
h All geometries optimized. 

c Calculations at SCF-optimized geometries ( S A / for biradicals and J A' for phenols). 
d Based on JANAF Tables and ref. 27. 


84 



Summary of previous results for C2H4 
• Energy map augmented by calculations of Dupuis et al. 
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